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We propose a new technique for comparing two Raman active samples. The 
method employs optical interference of the signals generated via coherent anti-Stokes 
Raman scattering (CARS) of broadband laser pulses with noisy spectra. It does not 
require spectrally resolved detection, and no prior knowledge about either the Ra- 
man spectrum of the samples, or the spectrum of the incident light is needed. We 
study the proposed method theoretically, and demonstrate it in a proof-of-principle 
experiment on Toluene and ortho-Xylene samples. 

I. INTRODUCTION. 

In the last decade, coherent nonlinear optical spectroscopy with femtosecond pulses has 
evolved into a powerful tool for chemical characterization, detection, and microscopy [1-8]. 
Ultrashort laser pulses exhibit high peak intensities, which result in high nonlinear signals, 
at low average power below the damage threshold of many systems of interest. Coher- 
ent anti-Stokes Raman scattering (CARS) with time-resolved (tr) and frequency-resolved 
(fr) detection has been recently complemented by a number of methods employing shaped 
femtosecond pulses. By adjusting the amplitudes and phases of the spectral components 
of a broadband pulse, one is able to enhance Raman excitation of one molecular species 
while not exciting another, thus gaining sensitivity to chemical structure [9] . Pulse shaping 
also proved useful in suppressing the non-resonance background signal which usually reduces 
CARS sensitivity [10-13]. Pulse shaping approaches to coherent nonlinear spectroscopy typ- 
ically rely on the availability of broadband pulses with smooth well-characterized spectral 
and temporal profiles [14]. 

In many situations, the goal of the spectroscopic analysis is a quick "yes" or "no" answer 
to the question of whether the two samples of interest are similar. One example is a quality 
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control task in which the sample in question is compared with a reference. Here, it is desirable 
to avoid experiments which require time-consuming scanning procedures, e.g. delay scanning 
in tr-CARS or frequency scanning in fr-CARS detection. Single-shot techniques based on 
femtosecond pulse shaping, although very quick, are often sensitive to the a priori knowledge 
about the anticipated spectral response from the sample of interest. Another quick scan- 
less approach to molecular analysis known as multiplex CARS (broadband excitation with 
narrowband probing) relies on spectrally resolved detection and analysis of the signal. 

In this paper we propose and demonstrate a technique for the direct spectroscopic com- 
parison of two Raman active samples. In the following text we refer to them as "reference" 
("i?") and "sample" ("aS"), and imply that the goal of the experiment is to establish the 
degree of their similarity. The proposed method does not require time or frequency scanning, 
and can be implemented without technically involved pulse shaping. Spectral analysis of 
the signal is not required; rather, an integral power of a single anti-Stokes beam is detected. 
The scheme does not rely on the a priori knowledge of Raman spectra of either the sample 
of interest or the reference. Moreover, it requires neither characterization nor control over 
the spectrum of the input beams. One can employ probe pulses with unknown random 
spectral profiles, as long as the latter exhibit sufficiently narrow features. One example of 
such spectrum is that of a transform-limited pulse sent through a random scatterer. 

II. METHOD 

The idea of the method is shown in Fig. 1(a). Laser pulses pass through the reference 
and sample - illustrated in the Figure by two glasses of wine, - generating a nonlinear 
spectroscopic response in both of them. If the two materials are similar (different), their 
spectroscopic responses are similar (different) as well. Varying the phase $ added to the 
nonlinear optical signal between the two media, one gains 100% interference contrast in the 
case of identical samples ("5 = i? "), and a lower contrast if S and R differ ("5 7^ i? "). As 
shown below, the interference contrast is a natural measure of the similarity between the 5' 
and R spectra. 

We implement the idea described above using femtosecond CARS spectroscopy. In fem- 
tosecond CARS, broadband pump and Stokes laser pulses excite a long-lived Raman coher- 
ence in the medium. A probe pulse, which can be either narrow- or broadband, stimulates 
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FIG. 1: Color online. The idea of the method, (a): Two nonlinear signals (thin red and yellow 
arrows) are generated in two media of interest, Sample and Reference, and interfere at the detector. 
Variable phase shift $ is added to one of the signals. Interference is strongest when the sample and 
the reference are similar, (b): Simplified energy-level diagram for the proposed method. Similar 
samples generate ansi-Stokes radiation (curly blue lines) at the same frequency, leading to strong 
interference. Chemically different samples generate signal at different frequencies, leading to no 
interference. 

radiation at the anti-Stokes frequency ujas = — + ^pr^ where cj^, ujs^ oJpr are the 
frequencies of the pump, Stokes, and probe spectra, respectively. The anti-Stokes signal 
consists of two parts. The resonant signal, whose time duration is determined by the life 
times of the Raman modes, is given by [15-17] 

1 



(1) 



Here n enumerates Raman resonances (e.g. vibrational energy levels of the molecule), 
are their energies, are their widths, 

/oo 
du'E*{uo' -n)E{uj') (2) 
-CXO 

is the two-photon excitation spectrum, and E{u)) is the spectral amplitude of the joint pump. 
Stokes, and probe field. The non-resonant background (NRB) signal, corresponding to the 
instantaneous electronic responce, can be approximated as 

/CXO 
Eicj - n)Ain) dn (3) 
-CXO 

and can be either weaker or stronger than the resonant one. In order to separate from 
Ejir^ one can delay probe pulses with respect to pump and Stokes pulses, as described below. 
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Consider, first, a model situation in which both R and S Raman spectra consist of 
a single line of width F, centered at and fi^, respectively (Fig. 1(b)). The idea of our 
method is most transparent in the multiplex CARS scheme, where a broadband pump-Stokes 
excitation is followed by a narrowband probe. For the latter, we assume that Epr{uj) = Eq 
within a narrow interval duj near the probe frequency cj^^, and Epr{uj) = otherwise. By our 
design, the CARS signal acquires an extra phase factor e^^ after exiting the sample medium. 
Assuming for simplicity that the broadband excitation is uniform, i.e. A{Q) — Aq in the 
relevant frequency range of the excited Raman modes, we obtain from Eq.(l) the following 
resonant signal generated in both materials: 

Cr 



Er{co) = 5UJ EqAq 



CO — LJr, 



(4) 



CARS signal from the sample is in the vicinity of uj^j^j^^ — Qs + ^pr while for the reference 
it is near ujcars — + ^pr- When the phase $ is scanned, these signals will interfere 
if fls — ^nd will not interfere otherwise (Fig. 2). To quantify this interference, we 
calculate the integral signal by integrating Er{uj)E^{uj) over all frequencies. Introducing the 
dimensionless frequency w — {uj — ujpr)/T^ we have from Eq.(4) 
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W — Ws + I W — Wji + I 

where Wr^s = ^r,s/^^ and ()^ stands for the integration over frequency. The integral in 
Eq.(5) is readily taken by contour integration, which yields 

\CsCr 
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where 



-^multiplex = 27r5a;^ l-^oAp 

^5 — ^R 

Wrs = ^ , 

tan A = Wrs 

(t>c = argfQCij] . 



(7) 
(8) 
(9) 
(10) 



Equation (6) is the key for understanding the proposed method. If S and R are similar, 
then Wrs ~ 0, and the interference contrast in the integral signal is 100%: at $mm = ^ 
the resonant CARS intensity drops to zero, as seen in Fig.2(a,c). If the two samples are 



5 




FIG. 2: Color online. (a,b): Two-dimensional interference maps for the intensity of multiplex 
CARS probing two similar (a) and different (b) samples. (c,d): CARS signals integrated over 
frequency. 

completely different, the corresponding CARS signals are at two different frequencies, and 
therefore do not interfere ^ this is the situation shown in Fig.2(b,d). Finally, if S and 
R lines partially overlap, the integral interference pattern is different from that of two 
similar samples in two ways. First, its contrast is lower. If the S and R signals are of the 
same strength, Cs = Cr^ then the contrast is reduced due the factor 1/(1 + ujj^^) in the 
square brackets of Eq.(6). From Eq.(8), one can see that wrs niay indeed be viewed as a 
natural measure of the mismatch between the two Raman spectra. Second, the minimum 
and maximum of the integral signal are offset from their values $rnm = ^ and 

^max = by the value A + This situation is illustrated in Figs. 2 and 3. 

We now show that having a single narrowband probe pulse (as in multiplex CARS) is 
not essential for the proposed interference detection. Let us illuminate the S and R media 
with transform-limited broadband pump and Stokes pulses, and a probe pulse that consists 
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of N narrowband components with an average spectral width 5uj. The spectral separation 
between these components is not important: one may employ a set of well isolated lines, as 
well as a single broadband line with a random phase modulation under its spectral envelope. 
The only relevant parameter is the spectral correlation length Su of the complex probe field. 
It is important, however, that the phases of probe's spectral components are not correlated 
with each other. 

The resonant ant i- Stokes response now has the form 

N 



Er{uo) — Sou EiA{uJ — UJpri) 



1=1 



+ 



R 



UJ — UJpri ~ + 



(11) 



This equation is similar to Eq.(4), except that the single term EqAq is replaced by a sum 
of N complex amplitudes Ei — \Ei \ exp[i0^] multiplied by the Raman excitation amplitudes 
A{uj — UJpri)- The intensity of the CARS signal Er{uj)E^{uj) thus contains A^^ terms. 



Er{uj)E;{uj) = 5uj^ ^\EiEj\e' 



e'^Cs 



A{UJ — UJpri)A''{uJ — UJprj) 



Cr 



UJ — UJpri ~ + UJ — UJpri ~ ^R + 
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Or - iT 



, . (12) 

^prj — UJ — UJprj — — 

Within a given realization of a random probe spectrum, the interference pattern in the (a;, $) 
plane can be quite complex. CARS signal consists of many spectral lines, reflecting both the 
S and R Raman resonances convolved with the probe spectrum. Two examples are shown 
in Fig. 3. If the S = the two terms in the square brackets of Eq.(ll) are similar, and 
the interference contrast at each frequency is high (Fig.3(a)). However, if the two materials 
are different, the cross terms in Eq.(12) contribute to the interference pattern in a random 
uncorrelated way. Depending on the spectral amplitudes \Ei\ and phase differences 0^ — 0^, 
the interference fringes at different frequencies lose their contrast and cease to occur at the 
same values of $ . As a result, the intensity map in the (cl;, $) plane shows no regular fringe 
structure, as demonstrated in Fig. 3(b). 

The apparent randomness of the CARS intensity map disappears once it is averaged over 
several noise realizations. Since the phases 0^ of different spectral components of the probe 
fleld are not correlated, the cross terms in Eq.(12) average out, while those with i — j do not. 
The resulting expression is similar to Eq.(5), but with a normalization factor accounting for 
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FIG. 3: Color online. (a,b): Simulated interference map for CARS with noisy probe with (a): 
two Toluene samples, and (b): Toluene and ortho-Xylene samples. (c,d): CARS interference with 
a noisy probe, averaged over noise realizations and integrated over frequency, (c): two Toluene 
samples; (d): Toluene and ortho-Xylene samples. 

having spectral lines in the probe. By repeating the steps which lead to Equation (6), 
we obtain, for the case of a noisy probe pulse, an expression of exactly the same form. The 
only difference is the normalization constant A^noise • 

N 

N^oise = 27r5uj'\Ao\'J2{\Epr,^\') (13) 

i=l 

where (..) denotes averaging over noise realizations, and Aq is the Raman excitation prob- 
ability, which is assumed to be uniform in the relevant frequency range. If all the spectral 
lines of the probe pulse are of the same intensity, the integrated CARS intensity is simply 
times that generated by a single spectral component. 

If broadband pump and Stokes pulses excite several vibrational resonances, one has to 
add in Eq.(6) a sum over all S and R lines. For different substances, the lines will, in general. 
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be randomly placed with respect to each other. The value wrs (Eq.8) will be positive for 
some pairs of lines and negative for others. After summation over Qs and the minima 
and maxima of the resulting interference pattern may not shift much from $rnm — ^ and 
^max = 0. Yet the contrast of the interference pattern will remain low if the S and R 
materials differ significantly. The contrast depends on the average value of wrs ~ a natural 
measure of the similarity of the two spectra. 

We numerically simulated an interference of two noisy CARS signals from Toluene (CrHg) 
and o-Xylene (CgHio), whose Raman spectra are well known. The excitation pulse parame- 
ters were similar to those used in our experiments and described in Sec. III. Fig.3(a,b) shows 
two calculated interference patterns. In plot (a), Toluene was used as both the Sample and 
Reference, whereas replacing the Sample by o-Xylene resulted in plot (b). Integration over 
frequency produced the results which are plotted in panels (c,d) of the Figure and demon- 
strate a striking difference in the interference contrast for the cases of Toluene-Toluene and 
Toluene-o- Xylene interference. Within a wide set of parameters, the simulations showed a 
100% fringe visibility for the interference of two identical substances, and only ^17% visi- 
bility for a Toluene-o- Xylene combination, given an equal strength of their resonant signals. 
The ^17% visibility of the Toluene-o- Xylene interference is due to the partial overlap of 
their Raman spectral lines, and thus gives the measure of their spectral similarity in accord 
with Eq.(6). 

Our numerical tests allowed us to make several observations. First, the integrated pic- 
ture is insensitive to the spectral shapes of pump. Stokes, and probe pulses as long as the 
correlation length of noise introduced to the probe spectrum is smaller than the average 
separation between the S and R Raman lines. Second, having many independent lines in 
the probe spectrum is advantageous as compared to having a few lines. In the case of many 
lines, not only the signal is stronger (see Eq.(13)), but one also needs to average the inter- 
ference map over fewer noise realizations. Finally, we noticed that using random pulses for 
the excitation^ as well as for probing, leads to qualitatively similar results, although requires 
averaging over large numbers of noise realizations, and gives higher interference contrasts 
for the ^^S ^ i?" case. Assuming non-equal signal strengths further reduced the interference 
visibility, enabling an easy distinction between the similar-samples and different-samples 
combinations. 

As expected and confirmed by our simulations, non-resonant background increases the 
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fringe visibility even if the interfering signals are generated by two different materials. In- 
deed, NRB from Sample is identical to that from Reference, and their interference produces 
high-contrast oscillations with $. Hence, as in all CARS schemes, the proposed method 
requires suppressing the non-resonant signal as much as possible. This task is difficult to 
accomplish in multiplex CARS, where the overlap of all three excitation pulses in time is 
large. Ideally, one needs a long probe pulse of the duration comparable to the lifetime of the 
Raman coherence, but with a short front edge delayed in time with respect to the pump- 
Stokes excitation pulses [10]. Such an optical field can be generated by propagating an 
ultrashort pulse through a randomly scattering medium: the sharp front edge corresponds 
to the ballistic part, followed by a long tail of randomly scattered photons [18]. 

In the final series of calculations, instead of introducing random phase modulation in the 
probe spectrum, we used the transmission spectrum of a random layered medium, calculated 
using the transfer matrix method [19]. Our simulations show that the integral interference 
curves of Fig. 3 are well reproduced as long as the widths of well separated un-correlated 
lines in the probe spectrum are smaller than the differences between the Raman lines. We 
refer the reader to Ref.[20] for a detailed discussion of the transmission fines of a random 
layered medium, and of their dependance on the system parameters. 

III. EXPERIMENTAL PROCEDURE AND RESULTS 

Experimental results have been obtained using the setup shown in Fig.4. It consisted 
of a laser system based on a femtosecond TiiSapphire oscillator (Synergy, Femtolasers) , a 
regenerative amplifier (Spitfire Pro, Spectra Physics) and an optical parametric amplifier 
(OPA) (Topas, Light Conversion). The amplifier generated 3 mJ, 35 fs pulses at the central 
wavelength of 800 nm and 1 kHz repetition rate. A portion of the 800 nm beam was 
coupled into a home-built spectral shaper based on a 640-element liquid crystal spatial light 
modulator (CRI, USA). The shaper produced probe pulses with a random spectral profile. 
The latter consisted of 25 randomly distributed Lorentzians of 1 nm width (full width at half 
maximum) and random relative phases, as required in the proposed approach. One example 
of such a random probe spectrum is shown in Fig. 5 (a). The corresponding temporal envelope 
of the pulse is plotted in panel (b) , which demonstrates another key feature of the used pulse 
shape - a sharp rising edge and a long random pulse train behind it. 
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Another part (1 mJ) of the 800 nm beam was used to pump an OPA, which produced 
Stokes and pump pulses at 1240 (signal) and 1125 nm (second harmonic of idler), respec- 
tively. These pulses were synchronized in time, whereas the front edge of probe pulses was 
delayed by 200 fs with respect to the overlapping pump-Stokes pairs. All three beams were 
collimated in a vertical plane (see insert in Fig. 4) and focused with a 25 cm focal distance 
silver mirror into a 200 fim optical path cuvette with 50 /xm thick walls serving as Reference. 

The beams were then collimated by the second 25 cm focal distance silver mirror. Probe 
and Stokes beams, together with the generated CARS beam passed through a fixed 50 /xm 




FIG. 4: Experimental setup. Femtosecond Ti: Sapphire laser system with OPA generates pump 
pulses at 1240 nm (signal) and Stokes pulses at 1125 nm (second harmonic of idler). Probe pulses 
at 800 nm were coupled into a spectral pulse shaper. The shaper produced probe pulses with a 
random spectral profile. All three beams were collimated in a vertical plane (insert) and focused 
into the first (Reference) cuvette. The beams were then colfimated. Probe, Stokes and CARS 
beams passed through a fixed coverslip glass. Pump pulses, on the other hand, passed trough a 
similar coverslip glass mounted on a rotation stage. All beams were finally focused into the second 
(Sample) cuvette. CARS signal was spatially separated and coupled into a spectrometer. 
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coverslip glass. Pump pulses, on the other hand, passed trough a similar 50 /xm coverslip 
glass mounted on a rotation stage. The latter provided a variable phase shift for producing 
the interference fringes as required by our method. All beams were finally focused into the 
second cuvette with another 25 cm focal distance silver mirror. CARS signal was spatially 
separated and coupled into a spectrometer (Model 2035, McPherson) operating with the 
spectral resolution of 0.5 nm and equipped with a cooled CCD camera (iDus, Andor). 
Exposure time was set to 0.5 seconds, and CARS spectrum was recorded as a function of the 
coverslip angle 6. The energy of all input beams was set at 3 /xJ per pulse. Reference cuvette 




FIG. 5: An example of the probe pulse used in the experiment, a) Spectral intensity (solid black) 
and phase (dotted red); b) Temporal profile. 
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was filled with Toluene, whereas Sample cuvette contained either Toluene or o- Xylene. 

Fig.6(a,b) shows the recorded CARS interference signal for the Toluene- Toluene and 
Toluene-o- Xylene combinations. The pattern was averaged over five realizations of the 
random spectrum of probe pulses. High visibility interference fringes are clearly seen in plot 
(a), confirming that the two materials - Reference and Sample, are very similar in their 
CARS response. Unfortunately, the parabohc shape of the fringes in the ($, A) plane did 
not allow us to integrate the signal over frequency while preserving the fringe contrast. We 
attribute the fringe curvature to chromatic dispersion, which causes an additional phase 
accumulation between the 5' and R media, and cannot be easily compensated. Similar 
pattern was observed in the interference pattern for the non-resonant signals, and reproduced 
in our numerical analysis. We note that in the absence of dispersion, our calculations predict 
that the 2D interference map for the case of ''S = i?" is a series of vertical strips. Hence no 
frequency-resolved detection will be necessary to detect the fringe contrast in the absence 
of dispersion. 

Even though integrating the two-dimensional experimental signal over a;, and hence elim- 
inating the need for a spectrometer, was not possible at this stage, we demonstrate that the 
fringe visibility at each wavelength is quite strong. This is seen in Fig. 6(c), which shows 
two cross sections of the full map taken at two arbitrary values of A (two horizontal lines 
between plots (a) and (b)). In contrast, both the 2D interference map and the ID cross sec- 
tions for the Toluene-o- Xylene combination, shown respectively in Figs. 6(b) and (d), exhibit 
irregular fringes with much lower visibility. This qualitative result of the proof-of-principle 
experiment confirms the feasibility of the proposed method. 

IV. SUMMARY AND DISCUSSION 

We have shown theoretically, and demonstrated experimentally, that the interference of 
two CARS signals obtained by scattering noisy light from two Raman active samples offers 
a good measure of their chemical similarity. Provided the concentrations are similar, the 
interference is determined by a single parameter wrs (Eq.(8)), which represents the degree 
of overlap between the Raman lines, averaged over all pairs of lines. Even though spectrally 
resolved detection was required at the present stage of the experimental development, we 
plan to eliminate this requirement in the future work. This could be done by either compen- 
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FIG. 6: Color online. Experimental results. Two-dimensional interference patterns for Tolune- 
Toluene (a) and Toluene-o-Xylene combinations (b). (c,d): Respective cross-sections of the 2D 
maps at two fixed wavelengths denoted by dashed blue and solid red lines, values. 

sating the residual chromatic dispersion, or by filtering the detected signal with a band-pass 
filter, which was eflFectively demonstrated in this work. 

For accommodating the samples of significantly different concentration, our method will 
have to include an additional calibration step. Using the non-resonant background, the R 
and S contributions are first equalized by attenuating one of the two CARS signals until the 
contrast of the non-resonant interference reaches 100%. After that, NRB is eliminated by 
delaying probe pulses and the resonant interference signal is detected. Further, excitation 
in the fingerprint region, rather than across the whole Raman active range, can be used for 
increasing the detection sensitivity. Finally, the required noisy probe can be generated by 
sending a broadband pulse through a randomly scattering medium as discussed in the text 
above. With these improvements, the proposed interference method may become a powerful 
tool for a quick preliminary test: if the interference contrast is above a certain threshold, a 
more accurate frequency- or time-resolved analysis is executed. 

By exploiting the interference of two noisy pulses for retrieving the information about 
the samples in question, our method employs a general principle known as Coherence- 
Observation-by-Interference-Noise ("COIN") [21-23]. In COIN, the interference of two noisy 
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signals is used to deduce the degree of their coherence. In our case, the latter is equivalent 
to the degree of their spectral similarity. If the two spectra of interest are similar, the two 
signals (however noisy!) will interfere and produce fringes of high contrast. When, on the 
other hand, the spectra are different, the interference is suppressed and the the fringe vis- 
ibility decreases. Random spectral noise has been recently exploited in coherent nonlinear 
spectroscopy [20, 24, 25], and this work extends its apphcations to interferometric CARS 
techniques. 

This work was supported by DTRA, CFI and NSERC. 
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